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Abstract 

The systematic study of ground-state properties of the three-dimensional 
Fahcov-Kimball model is performed by a well-controlled numerical method. 
The results obtained are used to categorize the ground-state configurations 
according to common features for weak, intermediate and strong interactions. 
It is shown that only a few configuration types form the basic structure of the 
phase diagram. In particular, the largest regions of stability correspond to 
phase segregated configurations, striped configurations and configurations in 
which electrons are distributed in diagonal planes with incomplete chessboard 
structure. Near half-filling, mixtures of two phases with complete and incom- 
plete chessboard structure are determined. The relevance of these results for 
a description of real material is discussed. 

PACS nrs.:75.10.Lp, 71.27.+a, 71.28.+d, 71.30.+h 
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1 Introduction 

The Falicov-Kimball model (FKM) has become, since its introduction |T] in 1969, 
one of the most popular examples of a system of interacting electrons with short- 
range interactions. It has been used in the literature to study a great variety of 
many-body effects in rare-earth compounds, of which metal-insulator transitions, 
mixed-valence phenomena, and charge- density waves are the most common exam- 
ples j2]. The model is based on the coexistence of two different types of electronic 
states in a given material: localized, highly correlated ionic-like states and extended, 
uncorrelated, Bloch-like states. It is generally accepted that the above mentioned 
cooperative phenomena result from a change in the occupation numbers of these 
electronic states, which remain themselves basically unchanged in their character. 
Taking into account only the intra-atomic Coulomb interaction between the two 
types of states, the Hamiltonian of the spinless FKM can be written as the sum of 
three terms: 

ij i i 

where f^, fi are the creation and annihilation operators for an electron in the local- 
ized state at lattice site i with binding energy Ej and df , di are the creation and 
annihilation operators of the itinerant spinless electrons in the d-hand Wannier state 
at site i. 

The first term of (1) is the kinetic energy corresponding to quantum-mechanical 
hopping of the itinerant d electrons between sites i and j. These intersite hopping 
transitions are described by the matrix elements tij, which are —t if i and j are the 
nearest neighbours and zero otherwise (in the following all parameters are measured 
in units of t). The second term represents the on-site Coulomb interaction between 
the ci-band electrons with density Ud = Nd/ L = j^J^i dfdi and the localized / elec- 
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trons with density Uf = Nf / L = \ Y^i ft ft: where L is the number of lattice sites. 
The third term stands for the locahzed / electrons whose sharp energy level is Ef. 

Since in this spinless version of the FKM without hybridization the /-electron 
occupation number ftfi of each site i commutes with the Hamiltonian (1), the /- 
electron occupation number is a good quantum number, taking only two values: 
Wj = 1 or 0, according to whether or not the site i is occupied by the localized / 
electron. 

Now the Hamiltonian (1) can be written as 

H = Y,hi,dtdj + EfJ2wi, (2) 

ij i 

where hij{w) = tij + Uwi6ij. 

Thus for a given /-electron configuration w = {wi,W2---Wl} defined on the 
three-dimensional lattice with periodic boundary conditions, the Hamiltonian (2) is 
the second-quantized version of the single-particle Hamiltonian h{w) = T + UW, so 
the investigation of the model (2) is reduced to the investigation of the spectrum of 
h for different configurations of / electrons. 

Despite its relative simplicity and an impressive research activity in the past, the 
properties of this model remained unclear for a long time. The crucial break in this 
direction has been done recently by exact analytical O IH El and numerical [HI U\ 
calculations. These calculations showed that the spinless FKM can describe (at least 
qualitatively) such important phenomena observed experimentally in some rare-earth 
and transition metal compounds like the discontinuous valence and metal insulator 
transitions, phase separation, charge ordering, stripes formation, etc. In addition, it 
was found jH] that at non-zero temperatures the model is able to provide the qualita- 
tive explanation for the anomalous large values of the specific heat coefficient and for 
the extremely large changes of the electrical conductivity found in some intermediate 
valence compounds (e.g., in SmBg). These results indicate that the spinless FKM, 
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in spite of its simplicity, could be a convenient microscopic model for a description 
of ground-state, thermodynamic and transport properties of real materials. How- 
ever, real materials are usually three dimensional while the most of above mentioned 
results have been obtained for the limiting cases of D = 1,D = 2 and D = oo. 
Thus one can ask if these results, or at least some of them hold also in three di- 
mensions. This is the question that we would like to answer in this paper. Here we 
focus our attention on the ground-state properties of model. The special attention 
is devoted to examine the three dimensional analogs of phase segregation, charge 
ordering, stripes formation and metal-insulator transitions observed in D = 1 and 
D = 2. From this point of view the paper represents the first attempt to describe 
systematically the ground-state properties of the FKM in three dimensions. To at- 
tain this goal we use a well-controlled numerical method that we have elaborated 
recently jH]. The method is based on the simple modification of the exact diagonal- 
ization method on finite clusters and consists of following steps, (i) Chose a trial 
configuration w = {wi, W2 ■ ■ ■ wl}- (ii) Having w, U and Ej fixed, find all eigenval- 
ues Afc of h{w) = T + UW . (iii) For a given A''^ = Wi determine the ground-state 
energy E{w) = Ylk=i^ + EjNj of a particular /-electron configuration w by fill- 
ing in the lowest Nd = L — Nf one-electron levels (here we consider only the case 
Nf + Nd = L, which is the point of the special interest for valence and metal-insulator 
transitions caused by promotion of electrons from localized / orbitals (/" /"~^) 
to the conduction band states), (iv) Generate a new configuration w' by moving 
a randomly chosen electron to a new position which is chosen also as random, (v) 
Calculate the ground-state energy E{w'). If E{w') < E{w) the new configuration 
is accepted, otherwise w' is rejected. Then the steps (ii)-(v) are repeated until the 
convergence (for given U and Ef) is reached. Of course, one can move instead of one 
electron (in step (iv)) simultaneously two or more electrons, thereby the convergence 
of method is improved. Indeed, tests that we have performed for a wide range of the 
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model parameters showed that the latter implementation of the method, in which 

1 < p < Pmax electrons {p should be chosen at random) are moved to new positions 
overcomes better the local minima of the ground state energy. In this paper we 
perform calculations with Pmax = Nf. The main advantage of this implementation is 
that in any iteration step the system has a chance to lower its energy (even if it is in a 
local minimum), thereby the problem of local minima is strongly reduced (in princi- 
ple, the method becomes exact if the number of iteration steps goes to infinity). On 
the other hand a disadvantage of this selection is that the method converges slower 
than for Pmax = 2 and Pmax = 3. To speed up the convergence of the method (for 
Pmax = Nf) and still to hold its advantage we generate instead the random number 
p (in step (iv)) the pseudo-random number p that probability of choosing decreases 
(according to the power law) with increasing p. Such a modification improves consid- 
erably the convergence of the method. Repeating this procedure for different values 
of Ef and U one can immediately study the dependence of the /-electron occupation 
number Nf = J2i w^T*™ '^^ ^^e /-level position Ef (valence transitions) or the phase 
diagram of the model in the nj — U plane. This method was first used in our recent 
paper jH] to study the ground-state properties of the one and two-dimensional FKM. 
It was found that for small and intermediate clusters, where the exact numerical 
solution is possible (L ~ 30), the method is able to reproduce exactly the ground 
states of the spinless FKM, even after relative small number of iterations (typically 
10000 per site). 

2 Results and discussion 

To examine ground-state properties of the spinless FKM in three dimensions we 
have performed an exhaustive numerical study of the model for weak {U = 1), 
intermediate {U = 2) and strong {U = 8) interactions. For each selected value of U 
and Nf [Nf = 0, 1, ..,L) the ground-state configuration y;™" is determined by the 
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above described method (we remember that the total filhng is fixed at 1). To reveal 
the finite-size effects numerical calculations were done on two different clusters of 
4x4x4 and 6x6x6 sites. A direct comparison of numerical results obtained 
on 4 X 4 X 4 and 6x6x6 clusters showed that the ground-state configurations fall 
into several different categories which stability regions are practically independent of 
L. Let us start a discussion of our results with a description of these configuration 
types for different values of U and Nf (in the remainder of the paper the values of 
Nf always correspond to 6 x 6 x 6 cluster). 

The largest number of configuration types is observed in the weak-coupling limit. 
Going with A'^^ from zero to half-filling {Nf = L/2) we have observed the following 
configuration types for U = 1. At low /-electron concentrations the ground-states 
are the phase segregated configurations (/-electrons clump together while remaining 
part of lattice is free of /-electrons) listed in Fig. la for two selected values of Nf. 
Since the ground-states corresponding to the segregated configurations are metal- 
lic we arrive at an important conclusion, and namely, that the metallic domain 
that exists in the one and two dimensional FKM persists also in three dimensions. 
In the one dimensional case the region of stability of this metallic domain was re- 
stricted to low /-electron concentrations Uf < 1/4 and small Coulomb interactions 
f/ < 1 jni • The numerical calculations performed in two dimensions revealed jH] 
that with increasing dimension the stability region of this metallic domain shifts to 
higher values of f/ {U ~ 3). From this point of view it is interesting to examine if this 
trend holds also for three dimensions. To verify this conjecture we have determined 
the ground-state configurations for increasing U at low /-electron concentrations on 
4x4x4, 6x6x6 and 8x8x8 clusters. We have found that the metallic region in 
D = 3 extends up to f/ ~ 5, what confirms the trend conjectured from two dimen- 
sional calculations (in addition, in accordance with two-dimensional results we have 
found that the critical value of Uc decreases with increasing Uf). It should be noted 
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that this result is crucial for description of insulator-metal transitions in real mate- 
rials (like rare-earth and transition metal compounds). In these materials the values 
of the interaction constant U are much larger than the values of hopping integrals 
tij PP, and thus for the correct description of valence and metal-insulator transitions 
in these compounds one has to take the limit U > t and not U < t. On the other 
hand it should be mentioned that in the Falicov-Kimball picture it is possible to get 
the metal-insulator transition much easier, for example by including spins. Indeed, 
numerical calculations performed for the spin-one-half FKM showed that the 
metallic domain is stable in this model for a wide range of model parameters, includ- 
ing large values of U and Uf. Above the region of phase segregation we have observed 
the region of stripes formation (A'^^ = 10, ..,20). In this region the /-electrons form 
the one-dimensional charge lines (stripes) that can be perpendicular or parallel (see 
Fig. lb). This result shows that the crucial mechanism leading to the stripes forma- 
tion in strongly correlated systems should be the competition between the kinetic 
and short-range Coulomb interaction. Going with Nj to higher values of Nf the 
stripes vanish and again appear at Nf = 26, however in a fully different distribution 
(see Fig. 2a). While at smaller values of Nf the stripes have been distributed inho- 
mogeneously (only over one half of lattice) the stripes in the region Nf = 26, .., 31 are 
distributed regularly. Above this region a new type of configurations (see Fig. 2b) 
starts to develop. We call them the diagonal charge planes with incomplete chess- 
board structure, since the /-electrons prefer to occupy the diagonal planes with slope 
1 and within these planes they form the chessboard structure. Of course, there is 
a considerable freedom in categorization of ground state configurations according to 
some common features and the case of diagonal planes used by us is only one of 
possible ways. The region of diagonal charge planes is relatively broad and extends 
up to Nf ~ 50. Then follows the region in which the chessboard structure starts to 
develop. As illustrated in Fig. 3a the /-electrons begin to occupy preferably the sites 
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of sublattice A, leaving the sublattice B free of /-electrons. In addition, the config- 
urations that can be considered as mixtures of previous configuration types are also 
observed in this region (see Fig. 3b). However, with increasing Nf the configurations 
of chessboard type become dominant. Analysing these configurations we have found 
that the transition to the purely chessboard configuration realizes through several 
steps. The first step, the formation of the chessboard structure has been illustrated 
in Fig. 3a. The second step is shown in Fig. 4a. It is seen that the chessboard struc- 
ture is fully developed in some regions (planes) that are separated by planes with 
incomplete developed chessboard structure. Such a type of distribution is replaced 
for larger values of A'^^ by a new type of distributions (step three), where both re- 
gions with complete and incomplete chessboard structure have the three-dimensional 
character (see Fig. 4b). 

The same picture we have observed also for intermediate values of Coulomb in- 
teractions {U = 2). The larger values of U only slightly modify the stability regions 
of some phases, but no new configuration types appear. In particular, the domain 
of phase segregation, as well as the domain of stripes formation are reduced while 
the domain of diagonal planes with chessboard structure increases. This trend is 
observed also for larger values of U. In the strong coupling limit {U = 8) the phase 
segregated and striped phases absent and the region of stability the diagonal planes 
extends to relatively small values of Nf ~ 20. Below this value a homogeneous distri- 
bution of /-electrons is observed. Thus we can conclude that all fundamental results 
found in one and two-dimensional solutions of the FKM (the phase segregation, the 
stripes formation, the phase separation, etc.) holds also in three dimensions, thereby 
the FKM becomes interesting for a description of ground-state properties (e.g., va- 
lence and metal-insulator transitions induced by doping and pressure) of real (three 
dimensional) systems The work in this direction is currently in progress. 

In summary, the ground-state properties of the three-dimensional FKM were 
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examined by a well-controlled numerical method. The results obtained were used to 
categorize the ground-state configurations according to common features for weak 
[U = 1), intermediate {U = 2) and strong interactions {U = 8). It was shown 
that only a few configuration types form the basic structure of the phase diagram 
in the nj — U plane. In particular, the largest regions of stability correspond to 
phase segregated configurations, striped configurations and configurations in which 
electrons are distributed in diagonal planes with incomplete chessboard structure. 
Near half-filling, mixtures of two phases with complete and incomplete chessboard 
structure were determined. 

This work was supported by the Slovak Grant Agency VEGA under grant No. 
2/4060/04 and the Science and Technology Assistance Agency under Grant APVT- 
20-021602. Numerical results were obtained using computational resources of the 
Computing Centre of the Slovak Academy of Sciences. 
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a) 




N.= 12 N.= 18 



Figure 1: Typical examples of phase segregated (a) and striped (b) configurations 
obtained for U = 1 and L = 6 x 6 x 6. Large dots: occupied sites; small dots: vacant 
sites. 
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Figure 2: Typical examples of striped configurations with regular distribution (a) 
and diagonal charge planes with an incomplete chessboard structure (b) obtained for 
U = 1 and L = 6 X 6 X 6. 
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Figure 3: The ground state configurations for intermediate /-electron concentar- 
tions. (a) Tlie formation of the chessboard structure, (b) The examples of ground- 
state configurations that can be considered as mixtures of configuration types with 
smaller nj ([/ = 1, L = 6x6x6). 
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Figure 4: Examples of an incomplete chessboard structure obtained for U = 1 and 
L = 6 X 6 X 6. (a) The chessboard structure is fully developed in some regions (planes) 
that are separated by planes with incomplete developed chessboard structure, (b) 
Both regions with complete and incomplete chessboard structure have the three- 
dimensional character. 
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